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In a pilot study, we use the topological charge density defined by the eigenmodes of the overlap 
Dirac operator (with ultraviolet filtering by mode-truncation) to search for lumps of topological 
charge in SU{2) pure gauge theory. Augmenting this search with periodic and antiperiodic temporal 
boundary conditions for the overlap fermions, we demonstrate that the lumps can be classified 
either as calorons or as separate caloron constituents (dyons). Inside the topological charge clusters 
the (smeared) Polyakov loop is found to show the typical profile characteristic for calorons and 
dyons. This investigation, motivated by recent caloron/dyon model studies, is performed at the 
deconfinement phase transition for SU{2) gluodynamics on 20'' x 6 lattices described by the tadpole 
improved Liischer-Weisz action. The transition point has been carefully located. As a necessary 
condition for the caloron/dyon detection capability, we check that the LW action, in contrast to 
the Wilson action, generates lattice ensembles, for which the overlap Dirac eigenvalue spectrum 
smoothly behaves under smearing and under the change of the boundary conditions. 
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I. INTRODUCTION 

The two current confinement scena,rios, the nionopo- 
le ^^'^ vortex mechanism [4, 5] of confinement 

in SU{N) gauge theory have become unified within the 
Z{N) vortex picture [1, 0, i, i, [13, [HI- Yet there ex- 
ists the old hope to connect confinement also with the 
topological structure as understood in terms of instan- 
tons calorons [ll, dyons merons (l8| 

and more generic objects |19j . all being carriers of Pon- 
tryagin charge. ^ The main reason for this desire is to 
bring confinement, on first sight a rather abstract prop- 
erty of pure (lattice) Yang-Mills theory, in closer relation 
to the physical origin of chiral symmetry breaking and to 
the continuum theory. 

During the last decade new selfdual solutions have en- 
tered the discussion. The aim is now to explain confine- 
ment in such a model via a detour through finite temper- 
ature, at < T < Tdep- The new solutions are the Kraan- 
van Baal-Lee-Lu [20, l2l|, [23, 123] calorons with a general 
asymptotic holonomy Voo 4- ^'^^ necessarily in the 

center of SU{N). For SU{2) the asymptotic holonomy 
can be roughly identified with the real- valued spatial av- 
erage of the Polyakov loop L = (l/V) J^x'^^'^i^)- 

Very recently, Diakonov and Petrov have worked out 
a model [l^ based on a gas of interacting caloron con- 
stituents, i.e. selfdual dyons, that offers already a com- 
plete picture of confinement at finite temperature. Al- 
though the presence of both selfdual dyons and antisclf- 



^ For brevity, in this paper we understand "calorons" , "dyons" 
and "selfdual" as including also "anticalorons" , "antidyon" and 
"antiselfdual" . 



dual antidyons has been ignored so far, the model is a 
convincing step forward. For the success of this descrip- 
tion (which describes also the limit of low temperature !) 
the maximally nontrivial holonomy of the gauge field is 
crucial. 

In one paper [l^ , authored last year at Humboldt Uni- 
versity, the capability of a caloron gas model to explain 
confinement has been explored in a Monte Carlo study. 
In this caloron model the opposite extreme case of dyons 
bound in calorons is dealt with. The importance of max- 
imally nontrivial holonomy for the correct choice of the 
caloron solutions to be used in the model was the central 
idea. Even a modest dissociation of calorons into slightly 
separated dyon-dyon pairs turned out sufficient to create 
a confining heavy-quark potential of the right order of 
magnitude. 

The assumptions of these models have to be confronted 
with the lattice. Will we ever have the chance to con- 
firm or disprove models of this type by analyzing generic 
Monte Carlo lattice configurations ? For some time al- 
ready our aim is to understand, albeit numerically, to 
what extent calorons and dyons coexist and are dis- 
cernible in the Euclidean gauge fields, most probably be- 
low T(jec- This has been the central question in our two 
previous lattice papers on the problem [l^, and is the 
central question also now. 

Traditionally, the presence of locally classical excita- 
tions like instantons in Monte Carlo lattice gauge fields 
has been explored by using methods like cooling [2^, [23, 
[H, [111, restricted cooling [30], smearing j3lJ, that re- 
placed the more demanding RG cycling [33, [11], and 
more general, by combinations of blocking and inverse 
blocking 34, 35, 36]. In the result, a well-localized topo- 
logical charge density (according to its field-theoretical 
definition [33, [11]) becomes visible. All these methods 
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actively change the gluonic field of the lattice configura- 
tions. Therefore they have been considered with some 
scepticism because of the methodical bias towards clas- 
sical fields. Concerning our particular point of view, one 
might criticize that these methods also tend to hide the 
presence of dyonic constituents inside calorons. 

For the evaluation of the topological charge density the 
situation has completely changed with the availability of 
methods dealing with overlap fermions 3^ 4^] - or other 
fermions with improved chiral properties 4l|, 42 1 - as a 



probe to explore gauge fields. A definition [43, |44[ of 
the topological charge density has been given involving 
the trace of the overlap Dirac operator. Then ultravi- 
olet filtering can be given a well-defined sense [3^, 
by restricting the trace to the lowest fermionic eigen- 
modes with |A| < Acut in the spectrum. The properties 
of this whole family of topological densities are strongly 
changing with Acut and have been investigated in detail 
in Ref. 

In our two lattice papers [U,!!^ on the caloron/dyon 
issue we were relying on the smearing technique and using 
the field-theoretical definition of the topological charge 
density in terms of an improved lattice field-strength ten- 
sor [48]. Additionally, in order to support the interpre- 
tation of the clusters of topological charge as calorons 
or dyons, the monopole content of the clusters has been 
analyzed in the maximal Abelian gauge. ^ 

With the present paper we return to the investigation 
of the caloron/dyon structure. We are replacing the field- 
theoretical topological charge density, which always re- 
quires smearing, by the overlap-based topological charge 
density in the ultraviolet filtered form mentioned above. 
We should remind the reader that chirally improved lat- 
tice fermions [4l|, \^ (another realization of Ginsparg- 
Wilson [sT] fermions) have already been used to ana- 
lyze unsmeared lattice configurations for the presence of 
calorons [s^]. Before that, unimproved Wilson fermions 
have been employed [53| for the description of nearly clas- 
sical calorons and dyons obtained by cooling. What was 
common to both techniques was inspired by the theo- 
retically known behavior of zero modes of caloron-like 
configurations Thus, particular emphasis was first 

given to the zero modes (or the real modes in case of the 
Wilson-Dirac operator), which must be present in con- 
figurations with topological charge Q 7^ 0, and to the 
effect on them of changin g th e boundary conditions for 
the Dirac operator [5^, 1531.1551 [56j . Confronting the zero- 
mode pattern with the picture revealed by smearing, it 
became clear [55| that the zero modes are part of the 
topological structure of a typical Monte Carlo configura- 
tion but cannot exhaustively explain it. 

In a recent paper |57l| reporting a collaborative project 
of the Humboldt University and Regensburg University 
lattice groups, it has been described how smearing and 



spectral filtering methods (with fermions and scalars) can 
be tuned to each other as far as the topological charge 
density is concerned. In our present context the relation 
between the fermionic filtering and the result of smear- 
ing is relevant. In the parameter space of competing 
methods (number of modes vs. smearing steps) a map- 
ping was defined by optimizing the cross-correlation be- 
tween the respective topological charge densities. As just 
two extreme examples we quote the observations that 10 
smearing steps are equivalent to the filtering by 50 modes, 
while 20 smearing steps are equivalent to not more than 
8 modes. These are only two arbitrarily taken cases of 
relatively mild and strong filtering. Of course, the struc- 
ture changes (the number of lumps decreases) with in- 
creasing smearing steps. Moreover, even the parameter 
mapping does not guarantee that the clusters of the re- 
spective densities exactly coincide. The pointwise overlap 
amounts only to 50 to 60 %. Although in Ref. [13] chi- 
rally improved fermions [U, [Z^] were employed instead of 
overlap fermions, the results give additional motivation 
and orientation for the present investigation and may be 
helpful to appreciate the new findings. We will explore 
the possibilities of identifying caloron-like and dyon-like 
structures for intermediate filtering employing 20 overlap 
eigenmodes. 

What is the conjectured physical picture ? Our previ- 
ous experience US EE] suggestSj_^ accordance with 



the model of Diakonov and Petrov [17| that a "plasma" 
including calorons (with nontrivial holonomy) and dis- 
solved dyonic constituents may describe the field struc- 
ture at T < Tdcc rather well. It fails, however, to de- 
scribe the essential features of lattice fields above Tdcc- It 
has been guessed that calorons with intermediate holon- 
omy would describe the topological structure in the 
high-T phase closely above Tdcc- A semiclassical evalu- 
ation of the path integral 60.] has shown, however, that 
the caloron becomes unstable against dissociation into 
dyons outside a narrow stability region \L\ > 0.72. On 
the lattice, by suitable measures of selfduality [13, [6l| , it 
has been observed that locally selfdual domains become 
suppressed above Tdcc- The topological susceptibility is 
known to slowly decrease (in the case of SU{2) gauge the- 
ory), and purely magnetic monopole excitations probably 
acquire an overwhelming importance. 

In the confined phase the caloron model indeed de- 
scribes confinement, even if the dissociation of calorons 
into dyons remains incomplete and within a description 
by a phenomenological choice of the p distribution. The 
size variable = d/(7rT) in the caloron case repre- 
sents a natural extension of the size parameter p^, usu- 
ally assigned to the (Euclidean) spherical lumps of ac- 
tion seen at T — > 0, to higher temperatures in the con- 



^ MAG was implemented on the lattice first in Refs. l49l |50(I . 



^ The usual Harrington-Shepard calorons [ssl , forming the basis of 
the first nonperturbative description of finite-T QCD [59l |. rep- 
resent the limiting case of trivial holonomy. 
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finement phase when the distance d between the con- 
stituents may be ^ 1/{ttT) or bigger. Dissociation is 
increasing with rising temperature towards T^cc- ^ This 
is in agreement with a practically temperature indepen- 
dent p-distribution for T < T^cc, similar to the usual 
instanton parametrization |15|]. Let us remark that at 
very low temperature the carriers of topological charge 
become more and more difficult to distinguish from gen- 
uine instantons by means of gauge-invariant observables 
alone (action and topological charge density). All this 
has lead to the conjecture that close to the deconfining 
phase transition the dyonic content of calorons becomes 
maximally manifest. Therefore we concentrate here first 
on this temperature. 

In the present investigation we have dispensed with (i) 
the use of smearing for the detection of clusters of topo- 
logical charge, and with (zi) the maximally Abelian gauge 
needed to determine the monopole content of the latter. 
We gave up, on the other hand, the exclusive focus on the 
zero mode(s) [13, of the configurations being under in- 
vestigation. We have concentrated instead on the effect 
of changing the fermionic boundary conditions on the 
whole overlap-based topological charge density mapped 
out by a given (not too large) number of eigenmodes. The 
dependence on boundary conditions has not yet been sys- 
tematically investigated. Our present paper is a first step 
in this direction and hopefully stimulates such a thorough 
investigation. Also concerning the 5*1/(2) gauge theory, 
this paper is the first application of the overlap-based 
topological charge density. For two colors it suffices to 
restrict oneself to the simultaneous consideration of peri- 
odic vs. antiperiodic boundary conditions. We shall see 
that the mode-truncated topological charge densities cor- 
responding to the two boundary conditions, completed 
by the local Polyakov loop, allow us to classify the vis- 
ible topological charge clusters as calorons and separate 
dyons, respectively. 

The paper is organized as follows. In section |ll] we 
define the technical details of our analysis: the ac- 
tion [H, , the overlap Dirac operator l39l.l40ll a nd the 
corresponding topological charge density 

~mm- Un- 
der conditions of confinement, but close to the transi- 
tion temperature, we critically check the stability of the 
fermionic topological charge (given by the index of the 
overlap Dirac operator) and the continuity of the low- 
lying spectrum under smearing and with respect to a 
change of fermionic boundary conditions. This check 
forces us to abandon the standard Wilson gauge action 
and motivates the choice of the Liischer- Weisz action that 
successfully passes the test. The bulk of investigations is 
performed using the tadpole-improved Liischer- Weisz ac- 
tion [gs'I. In section Hill we determine the critical /3imp,c 
for the quenched thermal phase transition with this ac- 
tion. The search for the transition is restricted to a 20^ x 6 



lattice that will be used in the following. Next, in sec- 
tion llVi we explain how the mode-truncated, overlap- 
based topological charge density obtained with the two 
different fermionic boundary conditions can be used to 
extract calorons and dyons from unsmeared configura- 
tions, but restricted to the resolution given by the num- 
ber of eigenmodes. This is practised for an ensemble 
generated on top of the phase transition. In the future 
we hope to proceed with this analysis deeper into the 
confinement and the deconfinement phases. A discussion 
of the results in the light of related work, our conclusions 
and an outlook will be presented in section IVl 



II. WILSON VS. LUSCHER- WEISZ ACTION: 
THE STABILITY OF THE DIRAC SPECTRUM 

A. The action 

We employ for the actual analysis of the caloron/dyon 
content of SU{2) gauge theory the tadpole improved ac- 
tion of the Liischer- Weisz form [6^, [63] 

S = Amp X! -^pi ~ X! ' 

pi rt 

where Spi and ^rt denote plaquette and 1x2 rectangular 
loop terms in the action, 

-^pi.rt = ^ Tr (1 - C/pi,,t) . (2) 

The parameter uo is the input tadpole improvement fac- 
tor taken here equal to the fourth root of the average 
plaquette Wiyi ^ ((1/2) Tr Upi). For SU{2) gauge the- 
ory, the tadpole factor uq has been selfconsistently deter- 
mined first in Ref. ^] for a few Amp values in the case 
of vanishing temperature on lattices with a suitable 
lattice size L for each value of the bare coupling constant. 
The result is given in Table I. For the convenience of the 
reader and later reference to the lattice scales we present 
the corresponding values of the string tension in lattice 
units. 

TABLE I: Details of the simulations with tadpole-improved 
Liischer- Weisz action at T = 



/3imp 


L 


Uo 


< P >^'^ 




2.7 


12 


0.87164 


0.87165(2) 


0.60(5) 


3.0 


12 


0.89485 


0.89478(2) 


0.366(8) 


3.1 


12 


0.90069 


0.90069(1) 


0.309(6) 


3.2 


16 


0.90578 


0.905765(3) 


0.258(5) 


3.3 


16 


0.91015 


0.910152(4) 


0.219(3) 


3.4 


20 


0.91402 


0.914020(2) 


0.180(3) 


3.5 


20 


0.91747 


0.917481(1) 


0.151(3) 



This tendency is also supported by the coohng results in Ref. |62|| . 
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In our simulations we have not included one-loop cor- 
rections to the coefficients nor considered non-planar 6- 
link loops the coefficient of which would be purely per- 
turbative. We have adopted the uq values obtained at 
zero temperature also for the simulations at T « Tjjec- 



B. The overlap Dirac operator 

The overlap Dirac operator is a particular solution of 
the Ginsparg- Wilson relation [FH 



Dj5 + -/5D^-Dj5D , 



(3) 



where p = 0(1) is a dimcnsionless parameter not to be 
confused ^ with the instanton or caloron size parameter 
above. These operators have nearly perfect chiral prop- 
erties. In particular, the Atiyah-Singer index theorem 
is fulfilled at finite lattice spacing, with N± clearly rec- 
ognizable zero modes with positive or negative chirality 
related to the topological charge 



Q 



index 



(4) 



This is unambiguous as long as the configurations satisfy 
certain weak smoothness requirements |68| . A Neuberger 
operator can be constructed starting from an arbitrary 
input Dirac operator (with bad chiral symmetry or with 
already improved chiral symmetry) through the steps we 
describe now. In our case, we take as the input kernel the 
simple Wilson-Dirac operator. In this case, the emerging 
Neuberger overlap operator [H, H^l is 



Do 



P 



1 



Dw/\f- 



Dw = M - 



(5) 

where Dw is the Wilson-Dirac operator with a negative 
mass term p/a. M is the Wilson hopping term with 
r ~ 1. An optimal choice is p « 1.4 . By construction, 
the operator £)ov satisfies the Ginsparg- Wilson relation. 
In order to compute the sign function in the alternative 
expression 



D 



D 



w = 75 sgn(i7w) , 



^^w = 75 D 



w 



6) 



we have used the minmax polynomial approximation [66 
Furthermore, the low-mode projection has been used: 80 
eigenmodes of the hermitean Wilson-Dirac operator TJ-yv 
have been treated explicitely. 

The topological charge density can be expressed in the 
form 43] 



q{x) 



-tr 



75 



2p 



(7) 



where tr denotes the trace only over color and spinor 
indices. This form of the topological charge density 
contains vacuum fluctuations of all scales. In the ap- 
parent chaos remarkable global, low-dimensional struc- 
tures [zOj [zHj [zl] are formed. They are three-dimensional 
at the percolation threshold |47l . The ultraviolet filtered 
(mode-truncated) density [45[ is written as a truncated 
sum over A as the dimcnsionless eigenvalues of a Dov/p, 



E 

|A|<A„t 



1 



A 



^AlsM^) , (8) 



however, shows clustering of topological c harg e |47| in 
four-dimensionally coherent clusters similar [57| to struc- 
tures usually revealed by smoothing the gauge field. This 
is the level of resolution where calorons and dyons may 
appear (or not). 

The integral over qx^^t i^) gives Q corresponding to the 
Atiyah-Singer index theorem, independent of Acut, be- 
cause only the zero modes contribute to Q according 
to their chirality. It is remarkable, but generally ob- 
served for the overlap Dirac operator, that if there are 
zero modes within a configuration, they all have the same 
chirality. 



C. The Dirac spectrum under smearing and 
varying boundary conditions 

The cross-relation between topological charge density 
and local structure of the Polyakov loop is typical for 
calorons and their constituents. In order to map out 
the Polyakov loop we need a modest amount of APE 
smearing in this study. A second role of smearing in 
our present context is that we want to monitor the in- 
dependence of the index of the overlap Dirac operator 
and a smooth dependence of the low lying spectrum on 
the number of APE smearing steps. We regard this as a 
necessary prerequisite that this part of the spectrum re- 
flects medium-scale and infrared properties only. Thus, a 
minimal requirement for the Dirac operator as well as for 
the lattice action (to prevent lattice artefacts that could 
give rise to unphysical zero modes) is the continuity of 
the spectrum under moderate smearing. This in fact se- 
lects admissible actions and an admissible range of the 
respective coupling. 

Smearing is an iterative sequence of four-dimensional 
link substitutions, where links are replaced by a weighted 
average of the links and the staples 



surrounding it: 



We prefer to keep this standard notation. 



(9) 

Here V denotes the projection onto the gauge group. 
For SU{2) this is just a rescaling of the matrix by a 
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scalar. We choose the smearing parameter a = 0.45 
foUowing |3h] where an optimal smearing schedule has 
been searched for. Within some limits, this parameter 
could be traded against the number of smearing steps. 
We allow for iVAPE < 10 iterations. 

It is known since long [t^ that the Wilson action is 
problematic with respect to dislocations. We remind the 
reader that the definition of a dislocation depends both 
on the action in use for the generation of gauge con- 
figurations and on the prescription chosen to define the 
topological charge (or charge density). In the case of 
the SU{2) Wilson action together with the (geometrical) 
Phillips-Stone topological charge Q , scaling of the topo- 
logical susceptibility [73] was an unsuspected fact before 
Pugh and Teper [7a| have shown that the observed value 
of (Q^) was dominated by excitations of size 0{a) that 
would not survive blocking. The lesson we draw from 
this example is that the bosonic topological charge 

Ogiuon = j d'^x ggiuon(a;) (10) 

to be assigned to a configuration by a suitably improved 
gluonic topological charge density 

'?gluon(2;) — 5'£'/ii/pcrTr {Ffii/Fpa) i (H) 

OZ TT 

and a geometrically defined version of the topological 
charge [zl, [z3| may be in systematical disagreement due 
to the presence of lattice artefacts. This can be the case 
even though it might not appear as a scaling violation 
of the topological susceptibility. In a first application of 
the overlap operator to SU{2) gauge fields generated with 
Wilson action [tI] , a reasonable continuum limit of the 
susceptibility {Q'^)/V with Q defined via the index of the 
overlap operator has been found. Therefore the Wilson 
action was not suspicious. The value of Xtopj however, 
was somewhat large compared to other estimates for the 
SU{2) gauge theory. ^ 

The field-theoretic definition of the topological charge 
density that we have used in our previous papers [2^ [25| 
employs the improved field strength tensor [48|. The ge- 
ometrical definition of the topological charge Q of a con- 
figuration is replaced in our present context by the index 
of the overlap Dirac operator, and the topological density 
by the corresponding expression ([7]) given above. At this 
point, checking the above-formulated requirements, dis- 
turbing features of the Wilson action are encountered. At 
first, the roughness of the configurations results in a rel- 
atively bad performance of the ARPACK package used 
to diagonalize the overlap Dirac operator. This proba- 
bly leads to a bad reproducibility of the measured index. 



The scaling property of the topological susceptibility would only 
be lost, resulting in a diverging susceptibility in the continuum 
limit, if local excitations would exist, that give rise to highly 
localized zero modes and would have a Wilson action less than 

sw < [zi. 



Thus, the latter can easily be misidentified due to zero 
modes pinned to dislocations. During the first smearing 
steps such dislocations become even more singular. 

As a typical example we show in Fig. [T] the lowest 20 
eigenvalues according to the two boundary conditions im- 
posed, without smearing and with 5 and 10 smearing 
steps, for a 20'^ x 6 configuration generated with the Wil- 
son action at /3 = 2.4 . This (3 was chosen below the 
critical value Pd^r = 6) = 2.4265(30) reported for the 
Wilson action in Ref. |80j- We see jumps of the measured 
topological charge, |A(5i„dex| = 1, between subsequent 
stages of smearing and occurring under a change of the 
boundary conditions (temporally periodic vs. antiperi- 
odic). During the first steps, smearing changes only the 
short range structure. The changing index counts here 
essentially the dislocations. Hence, the number of zero 
modes rapidly changes with the APE smearing steps. 

The fact that the Liischer-Weisz action is advantageous 
to facilitate our study, has been confirmed for a number 
of Amp values. The result is demonstrated in Fig. [2] for a 
typical configuration from a Liischer-Weisz ensemble at 
Amp = 3.2. The number of zero modes is independent of 
the type of temporal boundary condition and does not 
change with the number of APE smearing steps (as long 
as smearing is moderate, say iVAPE < 10). For Amp > 3.2 
we have never encountered such ambiguities as seen in the 
Wilson case. In Section IIIII we will see that the critical 
inverse gauge coupling for this action is Amp, c(A^t = 6) = 
3.248(2). The successful check presented in Fig. [5] has 
been performed for a situation close but clearly below 
the phase transition. 

We should stress, however, that configurations created 
by means of the Liischer-Weisz action may also turn out 
too "rough" at sufficiently low Amp values. For example, 
exploring the temperature range around Tdoc on a coarser 
lattice with Nt = 4 (i.e. at lower Amp), we found that 
the described ambiguities reappear. 

A surprising observation in the case of both actions 
is that the interval covered by the 20 lowest eigenval- 
ues does not systematically expand under the applica- 
tion of smearing steps. This differs from the behavior 
seen in Ref. [8l| for 16''' configurations generated with the 
(tree-level) Liischer-Weisz action. The spectrum there 
was considered not for the overlap Dirac operator but 
for the chirally improved Dirac operator proposed in 
Ref. [ill, H^l . It would be interesting to compare the two 
Dirac operators in their behavior under smearing for dif- 
ferent lattice ensembles (provided the smearing-induced 
changes are smooth). 



III. LOCATING THE FINITE TEMPERATURE 
PHASE TRANSITION 

The last observations make clear that we need to 
choose Nr > 6 for the purpose of this investigation. Let 
us now look for a more precise location of the deconfine- 
ment transition. On the 20'^ x 6 lattice, varying Amp, we 
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have studied the behavior of the Polyakov loop and of 
the Polyakov loop susceptibility. We used a polynomial 
fit for Mo as a function of /3imp, based on the measured 
values shown in Table I, in order to provide the corre- 
sponding tadpole improvement factor for each simulation 
point /?imp- We stress again that this non-perturbative 
determination, strictly speaking, is well-established only 
for temperature T = 0. 

We have measured the Polyakov loop and its suscep- 
tibility in the range from /3imp = 3.1 to 3.4 with differ- 
ent statistics per data point. The simulation data be- 
tween /?imp = 3.20 and 3.29, in the immediate vicinity 
of the phase transition, have been collected in 100,000 
to 300,000 Monte Carlo sweeps per Pi^p value while the 
Polyakov loop L = (l/V) Tr7'(x) was measured after 
every sweep. In the closer vicinity of the phase transi- 
tion we have fitted the susceptibility data by a Gaussian. 
The data and the fit of the susceptibility are presented in 
Fig. [3l For the determination of the errors, the blocked 
jackknife method was used with a block size of 2000 mea- 
surements. From the fit we are able to locate the de- 
confinement transition at /3imp,c = 3.248(2) for Nr — 6. 
This confirms our preliminary choice made in Sect. [Til of 
Amp — 3.2 for a check of smoothness of the overlap Dirac 
operator that should be done in the confinement phase on 
a lattice of the same size. Interpolating the data in Table 
I we estimate ^/aa = 0.236(5) at Amp c corresponding to 
rdcc/\/^ = 0.71(2). 



IV. FINDING CALORONS AND DYONS USING 
PERIODIC AND ANTIPERIODIC MODES 

In the Introduction we have argued why we should first 
search on top of the phase transition for calorons with 
nontrivial holonomy and why we anticipate to find them 
partly separated into dyons. We have chosen Amp — 3.25 
very close to the transition point for the following study 
of topological charge clustering. Our analysis is based on 
20 lowest-lying modes for an ensemble of 0(20) quenched 
configurations at the deconfinement transition. This was 
a realistic task within the capability of a standard mod- 
ern PC within a few weeks. 

The topological charge density of an equilibrium Monte 
Carlo field configuration is represented by the mode- 
truncated, i.e. ultraviolet filtered, topological charge 
density ([8]). In that definition the temporal boundary 
condition was not specified, that should be applied in 
the construction of the Wilson-Dirac and the Neuberger 
overlap o per ator From the work of Gattringer and 
Schaefer [54| we know that the single zero mode of a 
Q = ±1 Monte Carlo configuration eventually hops be- 
tween iVcoior positions. On the other hand, the topolog- 
ical charge density of a (quenched) lattice configuration 
cannot depend on the purely analysing fermions, in par- 
ticular not on the boundary conditions imposed to them. 
The most suggestive rule for the topological charge den- 
sity, if given by the zero-mode part of ([5]) alone, would 



be to average over the boundary conditions, which even- 
tually (but not always !) lead to a different localization 
of the zero mode. This recipe is now applied to the topo- 
logical charge density with the inclusion of the low-lying 
non-zero modes, too. 

We illustrate this in Fig. 2] for a classical charge Q = 2 
caloron solution with nontrivial holonomy in a state of 
maximal separation into four dyons. The upper pan- 
els show the gluonic definition (?giuon(2;) of the topolog- 
ical charge density and the profile of the Polyakov loop 
p{x) — (1/2) Ti'P{x) over a two-dimensional section of 
a 16"^ X 4 lattice. The gluonic topological charge den- 
sity recognizes all the four constituents as positive peaks 
while the Polyakov loop distinguishes the constituents ac- 
cording to the local holonomy, i.e. positive and negative 
values of the Polyakov loop. In the fermionic definition of 
the topological charge density, 9Aout(^)> content our- 
selves to only 20 lowest modes. We find that this filtered 
density depends on the boundary condition b, with b = p 
denoting periodic and b = a denoting antiperiodic tem- 
poral boundary conditions. The charge densities present 
a different profile depending on the type of boundary 
conditions. The antiperiodic boundary condition high- 
lights the constituents with negative local Polyakov loop, 
whereas the periodic boundary condition emphasizes the 
complementary constituents with positive local Polyakov 
loop. The "true" topological charge density (that is well- 
represented by the gluonic definition in this classical case) 
is well approximated by an average of the two fermionic 
topological charge density functions qx'\ (x) and <}^x \ (x), 

'^ilW-- E (i-y)<N5<'(-), (12) 

where the superscript b = p or b = a oi the modes (the 
subscript of the eigenvalues) refers to the boundary con- 
dition. 

Thus, for each boundary condition, we will search for 
peaks of the modulus of the corresponding fermionic 
topological charge density. In addition, in order to de- 
fine a size for the charge cloud surrounding the peaks, 
the respective topological charge density is separately 
subjected to a cluster analysis. As usual (see Ref. [23, 
m, 1131) the cluster analysis is a procedure to identify 
connected clusters among those lattice sites x G S, that 
have been selected by the condition that the modulus 
of the topological charge density \q{x)\ exceeds a certain 
threshold value gcut- Two sites x, i/ g 5, being neighbors 
on the lattice, belong to the same cluster, if the signs 
of q{x) and q{y) agree. Otherwise they belong to differ- 
ent clusters. Guided by Ref. [57|, the threshold is chosen 
relative to the maximal density in the configuration as 
9cut = I max2,(|g(x)|), safely above the point where the 
clusters coalesce and, finally, percolate. 

For the set C'^''-' {b = p, a) of clusters cf ^ in a configu- 
ration found by the cluster analysis of the two densities 
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^^'^ respectively we record the max- 

imal value of the modulus of the corresponding density, 



la^'^ I 

l^inax cluster, i I 



the sign sgn(g^^^ duster i) ^'^'^ corresponding space- 
time position X of the peak inside each cluster cf''^ . The 
main purpose of defining the finite size clusters around 
the peaks is to characterize the behavior of the Polyakov 
loop in the vicinity. The Polyakov loop is always mea- 
sured after A^ape — 10 smearing steps. Although the sign 
of the Polyakov loop at the cluster centers (topological 
density peaks) is found to be dictated by the temporal pe- 
riodicity/antiperiodicity imposed on the Dirac operator, 
the Polyakov loop is monitored all over the cluster to give 



= max iq 



(^)l, 



(13) 



auxiliary information. Its extremal values, P 



and 



-^min z' inside the clusters cf are recorded. At least one 
of the two corresponds to the fermionic boundary condi- 
tion that defines the clusters, being positive for the peri- 
odic boundary condition and negative for the antiperiodic 
boundary condition. 

Figs. [5] (a) and (b) show "cluster plots" representing 
two typical lattice configurations. A cluster cf'' of the 
topological charge density q'^^^x) is represented in the 
cluster plot by a filled circle {c^f^ ) for the periodic bound- 
ary condition or by a filled triangle {c^f'^ ) for the antiperi- 
odic boundary condition. The clusters c^^' are plotted in 
Figs. [5] at the appropriate position 



V^imax cluster, i 



cxtr cluster, i 



(14) 



in the (^max cluster ^ 



cxtr cluster 



plane. Here Pcxtr ciustor.i 
denotes either Pifi^ , or p£°), - , according to the p or a 
boundary condition that has defined the cluster through 
the corresponding topological charge density. Notice that 
this means that all circles appear in the upper and all 
triangles in the lower half-plane. 

Sometimes it happens that after changing the bound- 
ary condition from periodic to antiperiodic, one of the 

new clusters, cf'\ nearly coincides in its space-time posi- 

(p) 

tion X with one of the previous ones, , with a shift of 
the peak position less than a distance c? = 2a « 0.22 fm in 
space-time. This would correspond to the "not jumping" 
case of Ref. [s^ where, however, only the scalar density 
of a single zero mode was under consideration. In this 
case, such clusters, the circle c^^' and the triangle cf \ 
are connected in Fig. [5] by a broken line. The numbers 
close to the lines denote the approximate shift (0 or 1 or 
2) of the peak position. Such a pair represents a complete 
"caloron" , and the average over the respective topological 



charge densities q^^^x) and q^'^\x) locally represents the 
true topological charge density inside the caloron. For 
calorons the topological charge clusters defined for both 
types of boundary conditions are such that inside the 
clusters the extremal values of the Polyakov loop, P^^^ j 

and Pj^^nj, have clearly an opposite sign, indicative for 
the dipole structure of a caloron in terms of the Polyakov 
loop. 

Clusters that remained unpaired in this "cluster plot" 
have appeared only once, under only one type of bound- 
ary condition, such that the peak position could not be 
identified with a peak of the opposite boundary condi- 
tion, within a tolerance d < 2a. This corresponds to the 
"jumping" case of Ref. [13]. Such clusters do not have an 
obvious partner (with opposite sign Polyakov loop and 
same sign topological charge density) suitable to form a 
"caloron" . The length of the broken line attached to the 
unpaired filled symbols represents the difference between 
the maximum and the minimum of the Polyakov loop in- 
side the cluster. Numbers close to the unconnected lines 
denote the approximate distance (in the example, 3 or 4) 
between the cluster centers. In contrast to the caloron 
clusters both maximum and minimum of the Polyakov 
loop in an unpaired cluster are mostly of the same sign. 
In the few remaining cases the wrong-sign extremum is 
close to zero. Such clusters are called "dyons" because 
they, like the dyons in the classical Q — 2 caloron solu- 
tion shown in Fig. [31 are invisible to the fermions under 
the "wrong" boundary condition. 

The open circles around the filled symbols in the 
plots emphasize clusters which would have been localized 
knowing the zero mode(s) alone. These can also be clus- 
ters that we have to classify as calorons and as dyons. If 
they exist in the same configuration, these are clusters of 
a unique sign of the topological density, in accordance to 
the (yet unexplained) empirical fact that all zero modes 
of one configuration have the same sign of chirality. * 

Two characteristic objects that have been marked in 
Fig. [5] (b) as "CAL" and "D" are visualized in Fig. [6] 
in magnified form by their fermionic topological charge 
density profile q'^P^°-''{x) (left) and their Polyakov loop 
profile p{x) (right) within the occupied part of an x-y 
section: (a) for the caloron possessing the characteris- 
tic dipole structure of the Polyakov loop and (b) for the 
dyon possessing a broad maximum of the Polyakov loop. 
Let us stress that these objects have been identified in a 
generic Monte Carlo lattice configuration without cool- 
ing or smearing. To be sure, the Polyakov loop p{x) 
is presented after 10 smearing steps, which explains the 
smooth picture. 



For convenience we simplify the notation from now by dropping 
the subscript Acut from g^**' . 



* The cases of more than one zero mode per configuration were 
excluded from the analysis in Ref. [s^l . 
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We have also studied the relative separation of appro- 
priate dyon pairs. In Fig. [7] an histogram of dyon-dyon 
distances in our sample is presented. The first two bins 
correspond to calorons with distances d < 2a. The rest 
of the histogram with d > 2a contains pairs of suitably 
fitting dyon-dyon pairs, i.e. with the same sign of q{x) 
and an opposite sign of p{x), grouped into pairs accord- 
ing to the closest distance. The statistics does not war- 
rant so far the comparison with a specific model for the 
caloron/dyon plasma. 

Finally, we have collected in Fig. [8] all clusters of the 
whole sample of 20 configurations analogously to Fig. [5l 
The unpaired dyons are placed at their original positions 
in the (gmax duster, -Poxtr cluster) plane. The difference to 
Fig. [5] is that the two clusters close in space-time cor- 
responding to an undissociated caloron in this plot are 
re-located according to the average of the Polyakov loop 
assigned to the respective cluster as follows, 

cluster, i 2 y niax,i ' min,i J ' V-^^^ 

a value which is scattered around zero because of the 
dipole structure. Thus, in Fig. [51 each undissociated 
caloron is still represented by a close pair (with b — p 

and 6 = a) of open squares, now with |Pciustcr i\ < 0.25. 
The re-location according to the averaged Polyakov loop 
following Eq. (fT5|) leads in this scatter plot of clusters 
to a separation of points representing calorons and an- 
ticalorons (clustered in the ellipses) from the four types of 

dyons (clustered in the four circles) with |-Pciuster J > 0.5. 

The total number of isolated dyons (separated by a dis- 
tance d > 2a) is 113 plus 126 in this ensemble, whereas 
the number of dyons confined inside calorons (with a dis- 
tance d < 2a) amounts to 101 plus 101, meaning that on 
average approximately 10 calorons (dissociated or not) 
are present per configuration, if the resolution corre- 
sponds to 20 overlap eigenmodes. We emphasize again 
that the counting is a counting of peaks. All peaks get 
classified as dyons, regardless whether isolated or con- 
fined in calorons. The total number corresponds to a 
caloron (or dyon pair) density n^/^ — 230 MeV. This 
in the right ballpark set by the topological susceptibility, 
given the relative arbitrariness of the number of filtering 
modes. 



V. CONCLUSION 

In this paper we have continued our search for spe- 
cific KvB caloron-like features in finite-T lattice con- 
figurations. In a feasibility study we have for the first 
time employed overlap valence fermions for this diag- 
nostic purpose. More specifically, we have employed the 
dependence of eigenvectors and eigenvalues on the tem- 
poral boundary conditions imposed on the Dirac oper- 
ator that can be changed at will. Thereby, we have 
taken into consideration not only the zero mode(s) but 



the UV filtered topological charge density restricted to 
the 20 lowest modes per configuration. The dependence 
of the apparent caloron/dyon content on the number of 
eigenmodes has still to be systematically looked for. Ac- 
cording to Ref. [131 a resolution provided by 20 lowest 
fermionic eigenmodes, roughly corresponds to an amount 
of smoothing between 10 and 20 smearing steps. 

In our previous work (2^ . [25| we have used smearing 
and the corresponding gluonic topological density. The 
amount of smearing was, also somewhat arbitrarily, de- 
fined by the requirement that the string tension should 
not drop below 60 % of the full string tension [2^ which 
allowed for 50 or 25 ... 20 smearing steps in the confine- 
ment or deconfinement phase, respectively. Even more 
arbitrarily, the threshold for the definition of the clus- 
ters was set such that the density is split into a maximal 
number of clusters. Under these circumstances, a large 
number of shallow clusters entered the investigation be- 
fore only a small part of the clusters could be successfully 
characterized - by the monopole content - as calorons or 
dyons. 

In this work, apart from the number of modes dictated 
by the PC memory, we have fixed the cutoff (/cut in a re- 
gion where the number of clusters does not change with 
the cutoff and the size changes slowly. Moreover, the 
cluster centers were localized by the peaks of the mod- 
ulus of the fermionic topological density |9Aout(^)l ^'^d 
do not change anymore with the cutoff. Thus, the num- 
ber of clusters is determined essentially by the number 
of analysing modes that was adopted in anticipation of 
a physically acceptable density of dyon pairs. What we 
could show here is that with this resolution the cluster 
composition of the topological charge can be understood 
in terms of calorons and dyons without serious problems. 

All these clusters, once found, are seen to be accom- 
panied either by a dipole structure in the Polyakov loop 
p{x) or a broad maximum of the modulus of the Polyakov 
loop |p(a;)|. This shows that by means of the two topolog- 
ical densities (corresponding to periodic or antiperiodic 
temporal boundary conditions for overlap fermions) the 
task can be solved to identify calorons and dyonic con- 
stituents. 

In future investigations we will have to further spec- 
ify those conditions for filtering that make the cluster 
charges distributed around ±1 and ±1/2 , hopefully a 
very stable result. Furthermore, we hope for a better con- 
firmation of the caloron/dyon model by extending this 
study to lower temperature (where the model is good 
for describing confinement) and to study also the higher 
temperature region. 
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FIG. 1: The eigenvalue spectra of the Dirac overlap operator with (a) periodic and (b) antiperiodic temporal boundary 
conditions for equilibrium (A^ape ~ 0) and smeared (A^ape = 5 and 10) pure SU{2) gauge configurations on a 20'^ x 6 lattice, 
generated with the standard Wilson action at /3 = 2.40 . 
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FIG. 2: The eigenvalue spectra of the Dirac overlap operator with (a) periodic and (b) antiperiodic temporal boundary 
conditions for equilibrium (A'ape = 0) and smeared (A'ape = 5 and 10) pure SU{2) gauge configurations on a 20^ x 6 lattice, 
generated with the tadpole-improved Liischer-Weisz action at /3imp = 3.20 . 
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FIG. 3: (a) The ensemble average of the modulus \L\ of the average Polyakov loop L = y X)sTr7'(a;) as function of f3 for the 
tadpole-improved Liischer-Weisz action on a 20^ x 6 lattice, (b) The susceptibility of \L\ as function of f3. 
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(b) 

FIG. 4; (a) The gluonic topological charge density ggiuoii(a::) (left) and the Polyakov loop p{x) (right) for a classical Q = 2 
configuration generated at maximally nontrivial holonomy (asymptotically p{x) = 0) on a 16'^ * 4 lattice (with four dyons 
maximally separated in the (a;, 0)-plane). (b) The fermionic topological charge density q^^'''^\x) reconstructed out of the 20 
lowest eigenmodes of the overlap Dirac operator with periodic (right) and antiperiodic (left) temporal boundary conditions. 
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FIG. 5: The maxima of clusters of the fermionic |g(a;)i seen under periodic boundary condition (filled circles) and under antiperi- 
odic boundary condition (filled triangles) for two configurations (a) and (b) in the sample, shown in the (gmax cluster, -Pcxtr duster) 
plane (the precise meaning is explained in the text). Peaks at opposite-sign of Pcxtr duster, that are connected by dashed lines, 
have appeared under different boundary conditions at the same space-time position ("not jumping") and are interpreted as 
calorons. Isolated peaks have appeared only once under the respective boundary condition at the given position ("jumping") 
and are interpreted as dyons. The marked objects "D" and "CAL" in (b) are portrayed in detail in Fig. [B] 



16 




(b) 

FIG. 6; The fermionic topological charge density q^'''^°'\x) (left) and the Polyakov loop p{x) (right): (a) for a typical caloron 
cluster (when q'^'''\x) ~ q^°'\x)) and (b) for a typical dyon cluster (which was visible only in q''^\x)) from Fig. [5](b). The 
topological density and the Polyakov loop are represented as function over part of the (a;, y)-plane. Please notice the different 
scales for the topological charge density and for the Polyakov loop. The Polyakov loop is measured after A''ape ~ 10 smearing 
steps. 
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FIG. 7: The histogram of dyon-dyon distances in lattice units. The first two bins correspond to calorons which are unambigu- 
ously paired within distances d < 2a. The rest of the histogram with d > 2a refers to the remaining lumps grouped in suitable 
dyon-dyon pairs according to the closest distance. 
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FIG. 8: The topological clusters of the whole sample shown analogously to Fig. [5] The unpaired dyons are placed at their 
original (gmax cluster, fextr cluster) positious. The dyon pairs identified as calorons are finally re-located according to the average 
^'cluster over their original, opposite sign values Pmax and Pmin- 



